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1.  Introduction 

A  thin  spherical  layer  of  an  incompressible,  inviscid  fluid 
which  is  held  on  the  surface  S  of  a  rotating  ball  by  gravitation 
can  be  taken  for  some  purposes  as  an  approximation  to  the  Earth's 
atmosphere.   An  investigation  of  the  two  dimensional  vortical 
motion  in  such  a  layer  should  be  useful  for  the  understanding  of 
certain  observed  meteorological  phenomena.   For  example,  it  appears 
that  the  analysis  of  large  scale  closed  isobaric  systems  can  be 
based  on  a  knowledge  of  the  paths  of  concentrated  vortices. 

This  report  presents  a  study  of  the  vortical  motion  which  is 
due  to  the  existence  of  concentrated  vortices  (normal  to  S)  which 
are  confined  to  a  polar  cap  and  subject  to  a  boundary  condition 
along  a  circle  of  latitude. 

We  assume  that  the  departure  of  the  free  outer  surface  of  the 
layer  from  an  equilibrium  position  is  small;  and  that  the  tangen- 
tial acceleration  is  negligible  compared  with  the  Coriolis  force. 
We  also  assume  that  the  variation  of  the  Coriolis  force  with  lati- 
tude can  be  neglected.   In  other  words,  we  study  geostrophic  vor- 
tices on  a  sphere  as  contrasted  with  geostrophic  vortices  on  a 
plane.   The  latter  have  been  discussed  by  several  authors  in 
connection  with  the  motion  in  a  rotating  planar  tangential  layer 
as  an  approximation  to  the  motion  in  a  thin  layer  of  fluid  covering 
a  rotating  ball.   References  can  be  found  in  the  paper  by 
Morikawa,  [1]  , 

The  main  result  is  the  system  of  nonlinear  equations  for  the 
paths  of  the  vortices.  These  are  used  to  study  the  linear  stabil- 
ity of  the  motion  of  n  vortices  which  are  symmetrically  arranged 


along  a  circle  of  latitude  in  a  cap  with  or  without  a  polar  vortex. 
One  of  the  boundary  conditions  imposed  requires  the  velocity  to 
the  north  to  be  zero  along  the  circle  which  bounds  the  cap;  and 
the  other  requires  the  velocity  to  the  east  to  be  constant  along 
the  boundary.   The  investigation  is  restricted  to  analytical 
results  which  include  a  formula  for  linear  stability  for  the  case 
in  which  the  polar  vortex  is  either  held  fixed  or  else  has  zero 
strength.   Numerical  computations,  which  are  necessary  for 
decisions  about  the  linear  and  nonlinear  stability  of  given  con- 
figurations of  vortices,  will  be  discussed  in  another  report.   See 
Leiva,  [2]  .  „    .     , 


,1-l: 


2.  Equations  of  Motion 

Let  p  denote  the  distance  of  a  point  from  the  center  of  a 
ball  E  of  large  radius  a  which  rotates  with  constant  angular 
velocity  oj  about  a  polar  axis.   Let  (j)  ^'^'^  ^   denote  respectively 
the  longitude  and  the  colatitude  of  a  point  on  the  rotating 
spherical  surface  S  of  E.   Let  p  =  a  and  p  =  a+h((J),  0,t)  represent 
two  surfaces  which  contain  an  incompressible,  inviscid  fluid  which 
is  gravitationally  attracted  by  E.   Suppose  h  is  small  compared 
with  a. 

The  velocity  of  a  fluid  particle  relative  to  S  is  defined  by 
the  components 

u  =  (p  sin  6)   -5+  =  tangential  component  toward  the  east, 

d0 
V  =        -p  -rr-  =  tangential  component  toward  the  north. 


w  = 


-r£  =  radial  component. 


If  the  only  body  force  acting  is  that  due  to  the  gravitational 
potential  G  of  E,  then  the  basic  hydrodynamical  equations  which 
define  the  motion  of  the  fluid  relative  to  S  are  the  continuity 
equation 


^w  ,  „  ,   1 

P  -^  +  2W  + : 75- 

^   6d  sin   9 


and  the  momentum  equations 


^u  _    ^ ( V   s  in   6) 


=   0 


dU  UV     cot     e      ,      UW      ,       o  .  a  r^  a  1  1 

-__  _  +  —  +   2cow   Sin   6   -    2cDV   cos    6   =  -  -^ ^  .^    a  -r-i— 

dt  p  p  5    p   sin  0   ^^ 


dv  ,  vw  ,  u   cot  6  „^Q     1    1 

-jTT  +  —  +  +  2oou  cos  e  =  •= -^73—  ; 

dx    p       p  5  p  d0 

dw   V    u    „     .   a  1   ^1 

-r^    -    -  -  26DU  sin  e  =  -  •^—  -r . 

dT    p    p  5^  dp       - 

In  these  equations  the  differential  operator  with  respect  to  the 
time  means 

'^  =  ^  +    u    ^    _  y_    ^  +  w  ^ 
dTE"   "St   p  sin  0  "^"p'S^     '^  ' 

The  symbol  5   denotes  the  constant  density  of  the  fluid;  and  p 


2   2  2 

5    p  CD      sin   ^ 


stands   for  the   pressure    in 

P-L  =  P+5qG  - 

which  can  be  referred  to  as  a  modified  pressure. 

Since  h  is  supposed  to  be  small  compared  with  the  large 
radius  a  of  E;  and  since 

w((t),  e,  a,  t )  =  0 

let  us  neglect  the  radial  velocity  and  the  radial  variation  of  u 
and  V.   Let  us  assume  that  G  =  gp  where  g  is  a  constant;  and  that 

the  centrifugal  effects  manifested  by  the  partial  derivatives  of 

.       2   2      .    2o 
6    p  00      sm   0 

= can  be    ignored.      Let   us    also   assume   that   the  motion 

is  such  that  the  nonlinear  terras  in  the  tangential  momentum  equa- 
tions can  be  neglected;  and  that  the  radial  momentum  equation  can 
be    replaced  with   the   hydrostatic   law  .^^, 

Vi^,6,p,  t)    =   g6Q(h  +  a  -  p) 

which  satisfies  the  condition  that  the  pressure  is  zero  at  the  free 
surface  p  =  a+h. 


Under  the  above  assiimption  and  with  the  notation 


u((j),e,  t)  =  u(({),0,a,t)  , 


v((j),e,t)  =  v(({),e,a,t)  ; 


an  approximation  to  the  motion  is  determined  by  the  equations 


(2.1) 
(2.2) 

(2.3) 


^w  _|_     1 
"5p   a  sin  y 


du 


^(v  sin  6) 


^  -  2a)V  cos  0  = 


-W 


g    ^h 


0  ; 


asm 


■^  ^  ' 


^  +  2C0U  cos  0  =  I  ^ 


An  integration  of  (2.1)  gives 


w(^,0,a+h,t)  =^-^^ 


^(v  sin  6) 


SF 


The  kinematic  condition  at  the  free  surface  is 


w 


{^,d,  a+h,  t ) 


dh 


and  hence 

(2.4) 

dh 

dt   " 

h 

a  sin   0 

^(v  sin  0)    Su 
55       "^ 


We  suppose  now  that  h{^,6,t)    is  always  close  to  the  constant  value 


h  .   Then  if  we  introduce 


h  -  h 


Ti(<i),e,t)  =  — ^ 

a  linearization  of  (2.4)  yields 


o 


(2.5) 


1 

a  sin  6* 


^(v  sin  e)    ^u 


"5S       ^ 


while  the  momentum  equations  become 


(2.6) 
(2.7) 


du   _  ~ 

-^  -   2a)V  cos 


■^  +  2cDU  cos  e 


&K 


a  sin  0  ^  ' 


^  o  d 


a 
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Unless  stated  otherwise,  we  assume  in  what  follows  that  cu  is  not 


zero. 


3.  Geostrophic  Vortices 

If  we  neglect  the  variation  of  the  Corlolis  force  with  lati- 
tude and  take 


CO  cos  e  =  CO  cos  0^  =  O),  , 


then  the  last  equations  of  Section  2  reduce  to 


(3.1) 


(3.2) 


(3.3) 


'^t   a  sin  e 


du 

"5^ 


^v 

"3^ 


2a)-|^v  = 


o(v  sin  0)    ^u 

g^o   ^ 
a  sin  e  ^ 


gh^ 


-^-l-=-^|j 


Hereafter  we  confine  ourselves  to  a  study  of  these  equations.   As 
will  be  explained  In  the  sequel,  they  lead  to  what  are  called  geo- 
strophic vortices. 

The  elimination  of  t]  from  (3.2)  and  (3.3)  gives 


(3.4) 


h  1 

"St  ]  a  sin  e 


^(u  sin  6)        ^v 
5S       "5? 


=  20),  T] 


l''t  • 


The  quantity 


1 
a  sin  e 


^(u  sin  6)        hv 
3S "^  ^ 


Is  the  radial  component  of  vortlclty;  and  by  Integration  of  (3 .4) 
we  have 


(3.5) 


a  sin  e 


^(u  sin  6)        hv 
5^ ^  '^ 


2a)^Ti  +i:{^,9,0)  -  2cu-^Ti(({),0,O) 


This  implies  that  we  can  use  (3«5)>  (3  •2)  and  (3  •3)  as  a  basic  set 
of  equations  instead  of  (3.1);  (3 '2)  and  (3o3). 
Equations  (3.2)  and  (3.3)  imply 


(3.6) 

and 

(3.7) 


a^v  ^  u    2-   g^o 


^t' 


^""l  ■§?  "  sin  9   ^t^ 


■  2aj-^   ^     n2  ■ 
sin  0   f|  "^  "SOT 


If  these  are  used  to  eliminate  u  and  v  from  (3.5)^  we  find 


(3.8) 


sm 


O  r 

?r^  ^  ^^'^  ^  f^  +  — T7  u^  -  gh- 


sm  0  dtp      o 


^, .  1^ 


2o 
a  2CD-, 


;:Yr^  [(ii^,e,o)  -2a)^Ti((t),e,o)] 


p 
a  2a) 


gh 


i  C*(4>,e)  . 


o 


The  theory  of  Laplace  transforms  can  be  used  to  show  that  the 
steady  state  solution  of  (3.8)  is  such  that   ",   .'.;, 


(3.9) 


^t^oo'^^'^'^'^)  =^  ^(-f-.e) 


where   ip  must   satisfy 
1         ^ 


(3.10) 


sl^r^  -^e  ^'-''  "  "5^ 


M  +  _i_l!i  .  ^^^^1^ 


sm 


-e  h^' 


gh 


a-^i:*{^,6)    . 


The  function  ip   does  not  depend  on  t  and  the  associated  time  inde- 
pendent velocity  components  are 


8 


(J. 11)  u  =  L^^^u(|,e,t)  =i||; 

(3.12)  v  =  L,,^v(|.e.t)=^^^||. 

Note  that  (3<>10)  is  a  consequence  of  the  equations  (3.11),  (3  •12) 
and  the  vorticity  equation 


(^•^^)      a  sin  B 


h{U   sin  e)    ^V 
5^       ^ 


=  2cu,  L,  _Ti+C*((l),e) 


1  t^oo 

4a)-,  ^ 
^  o 

In  other  words,  (3.10)  is  implied  by  (3.5);  (3.2)  and  (3.3)  when 

we  ignore  inertial  forces. 

We  are  interested  in  the  solution  of  (3.10)  when  the 

vorticity  C^*{^,6)    is  constant  in  a  small  neighborhood  of  a  point 

{^*,6*)    and  zero  elsewhere.   This  solution,  as  is  well  known,  can 

be  derived  from  the  solution  of  the  idealized  equation  which  comes 

from  (3.10)  when  C,*i^,9)    is  assumed  to  represent  a  vortex  of 

strength  \i   concentrated  at  the  arbitrary  point  {^,6),      For  this 

reason  we  are  going  to  begin  with  the  mathematical  interpretation 

defined  by 

^2 
a'^[C('t'.e,0)  -2a)^Ti(c^,0,O)  =a.''Ci^,0)    =   -^—- 


aV5((|)-<t)-,)5(e-e.) 
a  Sin 


1 


where  5  symbolizes  the  Dirac  delta  function.   Hence  the  basic  equa- 
tion to  be  solved  is 


2  2 


D.14)  ^^^..ne^.—^^-^-.^-. s^^^ 


sm 


The  motion  in  a  thin  planar  layer  of  fluid  tangential  to  the 
surface  of  the  Earth  is  often  used  as  an  approximation  to  the 
actual  motion  of  the  Earth's  atmosphere  in  the  neighborhood  of  the 
point  of  tangency.   For  such  an  approximation  the  analogues  of 
(3.11);  (3.12)  and  (3.l4)  are  .  ' 


I  .' ', 


(3.15) 

(3.16) 
and 

(3.17) 


_1 
r 


SrU  ,  bV 


\J  =  h 


V 


1  ^X  . 


4cD-,x   |x5((|)-(l)-,  )5(r-r,  ) 

— vr-   + — 

gh  r- 


o 


1 


In  terms  of  the  polar  coordinates  (r,  (f)),  these  show  that  )(  must 
satisfy     ,,,.,_.,      ... 


The  only  physically  admissible  solution  of  (3.18)  is 


(3.19) 


M- 


y  =  -  T^  K 

^     2tt      o 


2C0-,  j—^ p 

- — ~  I  r  +  r^  -  2rr-,  cos  ((j)-i|)-,) 


/gh 


o 


where  K  [   1  denotes  the  zeroth  order  modified  Bessel  function  of 
o'- 

the  second  kind.   This  defines  what  is  called  a  geostrophic  vortex. 
The  motion  of  various  configurations  of  such  vortices  has  been 
studied  by  several  authors,  notably  G.  K.  Morikawa  [1],  [3]  whose 
papers  contain  detailed  explanations  and  references.   In  keeping 
with  what  seems  to  be  accepted  terminology  we  can  say  that  (3.1^) 


10 


defines  a  geostrophic  vortex  on  a  sphere. 

Our  object  now  is  to  investigate  the  nature  of  the  solution 
of  {j>,lk)    when  cd  is  not  zero.   This  equation  can  be  written  as 

sm  6  dtp  1 

where 

4a)-,  a 
(3.21)  v(v+l)  =  -  — ^  . 

^  o 

The  solution  of  (3.20)  can  be  inferred  from  the  solution  of 


This  equation  defines  the  Legendre  functions  of  degree  v.   Its 
general  solution  can  be  expressed  in  the  form 


f  =   c.P  (-cos  0)  + CgP  (cos  e) 


when  V  is  neither  zero  nor  an  integer  as  in  the  case  when  v  is 
defined  by  (3.21)  from  which 

V  =  -  -i  +  iq 
where 

i  =  y^  ; 


,  ,  iScjD-,  a 


and  CU-,  /   0.   In  the  neighborhood  of  the  north  pole  P  (cos  0)  is 
continuous  and  P  (l)  =  1;  but  in  this  neighborhood 


11 


^   (-COS  0)  =  P  -,    (-COS  e) 


behaves  like 


00 

2    ,     r       cos  qx  dx 

—  cosh  qw    / 

"^        ^Q  /2(cosh  X  -  cos  e) 


T.    r             n\        2  sm  VTT   „  r, 
P  (-cos  6)   ^  1^  .  i,n  6 

V  TT 


This  suggests  that  if  a  vortex  is  concentrated  at  {^-.,0^);    if  d, 
is  the  geodesic  distance  from  this  point  to  an  arbitrary  point 
((f),©)  on  S;  and  if  a-,  is  the  angle  a-,  =  d^/a,  then  the  fundamental 
solution  of  (3.20)  is 


•.. ) 


M- 


Sm  VTT 


P^(-cos  a^)  =  -^- 


Sm  TTV   V 


-COS  e  cos  6. 


-sin  0  sin  0-,  cos  (({)-(})-,) 


It  can  now  be  verified  by  direct  substitution  that 


(3.22) 


M- 


tl/    =   -r 7^ P  (-COS  a-.   ) 

^   4  sm  VTT  V       1 


does  indeed  satisfy  (3.20).   Furthermore,  computations  with  the 
velocity  components  f 


U 


tL 

l-a  sin  V 


—  -^  P   (-cos   a-,  )    ; 

TT    00       V  1 


V 


71 : 3 = ^TT    P     (-cos     O -,   )     ; 

4a  sm   0    sm  vtt   ^(|)     v  1'    ' 


show  that  in  the  neighborhood  of  {^^,6^),     (3 .22)  defines  a  vortical 
motion;  but  the  vortex  point  itself  remains  at  rest  —  it  possesses 
no  autonomous  motion.   In  fact  it  caji  be  shown  that 


12 


L  TT J=- P     (-COS     0-.   ) 

a-*oo   4    sm  vtt     v  1 

2CD 


li. 


K 


^      o 


/gh 


jr    +  r-,  -  2rr-,    cos    (({)-:|)-,) 


o 


Hence   we   can   say   that    (3.22)    represents   a  geostrophic  vortex  on   a 
rotating   sphere    in  the    same  way   as  we    say   that 


M- 


X  Ifjf 


K 


2co 


/ih 


i    /r^  +  r^  -  2rr      cos    ((ti-ct)^) 


represents  a  geostrophic  vortex  on  a  rotating  plane. 

For  an  arbitrary  distribution  of  n+1  vortices  on  the  sphere 
the  function  if/,    a  stream  function  is 


n 


1 

■d/  =  T, : T    M-  -P  (-cos  a  .) 

^       4  sm  VTT  4y^  J  V  ^      J ' 


where  a.  =  d ./a  and  d.  is  the  geodesic  distance  on  S  from  the 
point  of  concentration  of  the  j-th  vortex,  {^.,9    ),    to  an  arbitrary 

J     J 

point  ((f),  0).   That  is. 


cos  a-   =   cos  6   cos  0  .  +  sin  6    sin  9.    cos  ((])-({).) 
J  J  J         J 


The  associated  velocity  field  is  given  by 


U 


=  1  M  ^    1 4[y~:  i^-P   (-cos   o.)} 

a  d0        4a  sin  vtt   69  [  4^  ^j   v  ^  j  '  J 


V  = 


^ 


a 


sin  6  ^       4a   sin   9    sin  vtt   ^(J)  I  "tj^     J   v 


A{^   P^(.eo.  aj)} 


The  path  of  the  k-th  vortex  is  along  a  stream  line. 
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(^(x,y)  =  const),  and  its  velocity  is  equal  to  the  velocity  at 
(({'v'^if)  *^^^  ^°  ^-^^   °-^  ^^^  other  vortices.   Therefore  since 


U  =  a  sin  6 


^^ 


V  =  -a 


ae 

d^ 


the  equations  for  the  motion  of  i^-i^t^-i^)    are 


k'  k' 


and 


— ^ 

4a  sin  6    sin  vtt 


n 


4a  sin  0  sin  vir 


(f)=(f) 


k 


0=e,, 


If  the  motion  is  confined  to  a  part  D  of  the  sphere  S, 

boundary  conditions  must  be  imposed;  and  it  is  not  then  possible, 

in  general,  to  express  ?//   as  a  finite  sum  of  terms  of  the  type 

■n r-^^ P  (-cos  a.).   When  D  is  a  cap  bounded  by  a  circle  of 

4  sin  VTT  V  ^     .  J ' 

latitude  we  have  a  case  which  is  important  for  applications  to 
certain  meteorological  phenomena.   The  assumption  that  the 
Coriolis  force  is  independent  of  latitude  decreases  in  severity  as 
the  latitude  of  the  boimdary  circle  of  the  cap  increases. 


14 


k.    Geostrophic  Vortices  in  a  Cap 

Suppose  that  n+1  concentrated  geostrophic  vortices  exist  in 
the  cap  which  covers  the  north  pole  and  is  boi;inded  by  the  circle 
of  colatltude  0  =  A.   Suppose  that  along  this  boimdary  circle  the 
velocity  to  the  north  is  zero.   In  order  to  study  the  motion  of 
the  vortices  we  need  the  stream  function  which  satisfies 


(4.1) 


for 


_^^(sine)||(^,e)  . 


1   ^'^^ 


2 ^  +  V  {v+l)il/ 


=  1 


sin  9   b^ 
n  M.  .5  ((()-({)  .)5(e-e  .) 


J=0 


sin  8 , 


0   <   e,e  .    <   -K    ;         0   <   ^,^  .  ^  2v   ; 

J  J 


with  the  "boundary  conditions 


(4.2) 


V(<1),A)  =  i  ^^((t),A)  =  0 


or 


(4.3) 


^({j),A)  =  c  . 


The  solution  of  (4.1)  subject  to  (4.3)  can  be  expressed  in 
the  form 


1  -i^ 

(4.4)      ii/((t),e)  =bP    (cos  e) +  -;^-^-± >         ^L.P 


-COS  y  COS  y  . 

J 


-sin  0  Sin  0  .  COS  ((|)-(f)    ) 
J  J 


-,  n  oo 

+  -r — ^ \        \i  .  )         C    (m,A)F  (cos  e)F  (cos  6  .)cos  m((j)-(j)  .) 


j=0     ^   m=0 
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where  I^(x)  is  the  spherical  harroonic  of  degree  v  and  order  m. 
This  function  satisfies  ,  ,    ■  ,  i   ,  '   ,  ,   c 


dx  ^     ^  dx  V  ^  ' 


v(v+l)  - 


m 


1-x 


I^(x)  =  0 

V 


1    '. .  1, 


and  is  such  that 


V  ^  ' 


1  ,    m  =  0 


0  ,    m  7^  0  .   , 


p'^(x)  can  also  be  expressed  in  terms  of  P  (x).   We  have 
V  ^  V 


(4.5) 

The  expansion 


rs  ^  /o   d   P   (x) 

i^(x)  =  (-ird-x^r/^— ill . 

^  dx 


where 


and 


P   [-cos   Q   cos   e  .  -  sin  Q    sin  6  .    cos    (({)-(}>.)] 
V  J  J  J 

=  \         E    p"(cos   0  .)P^(-cos   0)    cos  m((})-(j).) 


Q)<e.<Q<ir;         (j),  i.    real 


^o   =   ^^ 


'm 


_    .    ^  ,in  2r(v  -m+1) 
"    ^~^>        r(v  +m+l)    ^ 


m   >  0  , 


P°(x)   =  p   (x)    ; 

V  V 


shows  that  the  boundary  condition  (4.3)  is  satisfied  if  we  take 
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h  =       C 

P  (COS  Aj 
V 

and 

P^(-cos  A) 

C  (rn.A)  =  -e  — 

•■  -"   P^(cos  A) 


m 


It  can  be  shovm  that  the  solution  (4.4)  is  unique;  and  that  it  is 
analytic  in  the  cap  minus  the  vortex  points. 

If  instead  of  (4.3)  we  impose  the  boundary  condition 

(4.6)  ^e(<l>,A)  =  c  , 

which  means  that  the  velocity  to  the  east  is  kept  constant  along 
the  boundary  circle  of  the  cap,  then  (4.4)  again  provides  the  solu- 
tion if  we  now  take 


and 


b   = 

=     - 

c 

sin 

AP^C 

cos   a) 

pB- 

(-cos   a) 

c 

V 

,(m. 

,A) 

V 
V 

(cos   a) 

where  p"^  (x)  denotes  the  derivative  of  E^(x). 

The  series  which  appear  in  (4.4)  can  be  summed  if  the  cap  is 
the  northern  hemisphere.   For  this  case  A  =  7r/2  and 

Cv^^'?)  =  +^m  • 
Then  since 

oo 

(4.7)     y^  e  F^(cos  e)P™(cos  0.)  cos  m(h-h.) 

=  P  [cos  6   cos  e  .  -  sin  6    sin  6.    cos  (i-i.)] 
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0   <   e,  e      <  TT   ;         e  +e      <  w   ;        ^,^ .    real   ; 
J  J  J 


we   see   that 


n 


(4.8)      f{i),e)=ToP   (cos0)  +  -r — r^ ^        ii.P 

^         /      Y'\T>     /  yv  /      4  Sin  VTT    •4— j^     J   V 


-cos  9  cos  6 


-sin  6  sin  9  .  cos{^-^ .) 

J  J 


1  " 

^    4     sin    VTT    4— rr    ^J     V 


J=0 


'   '  COS  9  COS  9  . 
3 


-Sin  0  sin  0  .  COS  (<!)-({).) 
3  <J    - 


where   for   the  boundary   condition 


^(i,|)   =  c 


the  upper  sign  is  to  be  taken  with 


P  (0)  ' 
V  ^  ' 


while  for 


fgi^,^)   =   c 


the  lower  sign  is  to  be  taken  with 


b  =  - 


Pv(0) 


Notice   that   P   [-cos  ©cos  6.-  sin0sin0.  cos  (({)-(]).)]    is    singular  at 

"  J  J  cl 

((f),  0 )  =  (|)  .,  e  . )  and  P  [cos0cose.-sinesin0.cos((|)-|).)]  is  singu- 

J     J  V  J  J  J 

lar  at  {^,9)   =  ^.,Tr-9.).      It  follows  from  the  form  of  (4.8)  that 
that  solution  can  be  constructed  from  the  fundamental  solution 


W^ 


-. P  [-cos  9   cos  9  .   -  sin  9    sin  9      cos  ((t)-<l)  •)] 


m  VTT   V 
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by  using  a  process  of  reflection  across  A  =  it/2.,  A  reflection 
technique  for  the  construction  of  the  solution  for  a  cap  other 
than  a  hemisphere  is  unknown  when  v   f^  0. 

If  V  =  0,  that  is  if  0)  =  0,  then  (4.1)  reduces  to 


(4.9) 


1    ^ 

sin  e  ^ 


,    ^2,     n  M.  .5((f)-(t).)5(e-e  .) 
.  2.  TJ^  -  ]—.  sin  B, 


This  equation  possesses  the  fundamental  solution 


J^  in   a? 


e 


e . 


tan     -^  +   tan     -^  -   2   tan  -^  tan  -^  cos    ( 


*-*j'] 


=  —-  Re    in  a 
2iy 


1   .^^ 
2 


e .    i(t) ; 


tan  -^  e""^-  tan  -^  ©   "^ 


Let  us  also  remark  that  if  we  employ  the  transformation 


/   4-    ^N  i't' 
z  =  (a  tan  -^je  ^  =  x  +  ly 

then  the  problem  of  solving  (4.9)  subject  to  0  _<  S  <  A  and  a 
boundary  condition  can  be  transformed  into  the  problem  of  solving 

l!i  +  l!|  =  ^il  6(x-x  )6(y-y.) 
^x^   ^y^   j=0  J     J      J 


/  P   2' 
subject  toO</x+y   <R  and  a  corresponding  boundary  condition. 

See  Peters,  [4] .   The  latter  equation,  as  is  well  known,  governs 

the  motion  of  rectilinear  vortices;  and  for  many  cases  it  can  be 

solved  by  using  reflection  techniques. 

The  velocity  field  defined  by  the  stream  function  (4.4)  is 

given  by 
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U  =   sin   0    (f)   =  - 


1   ^f 
a  "50 


4a     sin  vTT'Sin  6*^  = 


4    sin    VTT    -^ 


=  -4b 


sm  VTT'sm 


••p'  (cos   e) 


n 


+ 


4 — ;r       TV 

j=0       ^ 


-COS  0  COS  6  . 
3 


-sin  0  sin  6  .  cos  (({)-([) . ) 


sin  0  cos  0 


-cos  0  sin  0  .  cos  (<})-({).) 

J  J 


-    sm 


n  00  ,  ^ 

)        H.  EI  C„("^'^)P^    (c°s   e)F  (cos   0.)cos  in(^-(().) 


and 


V  =  -a  0   = 


1  ^ 

a  sin  0   ^ 


2  *         I  ^ 

-4a     sin  VTr-sin  0*0   =  4   sin  vtt  -^ 


n 


J^ 


H.P' 
J    V 


-cos    0    cos    0 . 

J 


-sin  0   sin  0.    cos    ((})-|).) 
J  J  . 


•  sin  0   sin  0.    sin    ((l)-ii).) 
J  J 


n  00 


-  y^  ^-  HI  ^    (m,A)P^(cos   0)P^(cos   0.)m   sin  m((t)-(|).)    . 


■po      ^   in=0 


The   following  equations    determine   the  motion  of   the   vortex 
at    {^]^>Sy.)' 


20 


(4.10)      4a     sin  vifsin  0    'L     = 


k   ^k 


-4b    sin  VTT   sin  6,  P   (cos   6,  ) 

k  V  k' 


n 


+  -^ —  ^.P' 

^T TT  1       V 


-cos  e,   cos  9  . 


-sin  e^  sine    cos  (<l'j^-(t'j) 


sin  6.   cos  e  . 


-cos  e^  sin  e     cos  ('j'j^-'t' J 


n  00  , 

-    sin  0     YZ  ^i  TZ  ^v^"''^)^    (^°s   ®k^^^'^°^    ^1^  cosin(|)^-(t)    )    , 


^   J^     J   fe     ^ 


^k   ^j 


■J  ■ 

(4.11)      -4a     sin  V7r*sin  S^^-e 


k     k 


n 


+  >         M-  -P 


-cos   0,    COS   6  . 


-sine.   Sin  e  .  cos  ((j)   -i  . ) 


.  sin  e,     sin  6.    sin    ((t)^-'!'-!) 
1^  3  ^      J 


-  yZ^iZZ  C^(m,A)F^(cos    ej^)l^(cos   6    )m  sinm(([)j^-(j)j )    . 


j^O     ^   m=0 


Suppose   now  that   one   vortex  of   strength  \i      is   at   the   north 
pole    and   that  n    (n  ^  2)    others   each   of  strength  \i,    are    symmetric- 
ally  situated   on   a  circle   of   latitude   within  the   cap.      In  other 
words,    suppose   that 


(^^,e^)  .  (^^,0)  , 


while 


%  =  ^'o 


(ij.e.).  ((j.i)^,^)^ 


j   =   1,  2,  . . . ,  n 


^J    =^ 


J 
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where  '  , . 

0    <  7    <   A    . 


For   this   configuration  the  velocity   to   the   north   of  the   polar 
vortex   is   determined  by  .       •      •  :     ■.  ' 


2  •  '^         t 

-4a     sin  V7r*0      =  \i  )         P   (-cos  -y )    sin  7«sin    ((f)   -^ .) 


n        CO  P^(cos   0,  ) 

V  ^  k^ 


-   11    L^  ^^  Y^  y^  C    (m,A)   -^:! : ^  P^(cos    0  .  )ra   sin    ( tb   -i.) 


Using    (4.5)   we   find 


2  •  n        , 

-4a     sin  vtt'6     =  \i  )         P   (cos  -y)    sin  7    sin    (({)  -(}).) 

J 

n 


-   ^^  ^  C^(l,A)P^(l)P^(cos  7)    sin  7    sin    i^Q-<^^) 


From  this  we  see  that 


e   =  0 
o 

.   .-   ■■.,■,  '•■■■■■   Ir  -'  -    -    ■        ■ 

because 


5~  sin  (4)^-i.)  =  >^  sin  [^^  -  (j-1)  ^]  =  0  . 

The  velocity  to  the  north  of  any  one  of  the  vortices  on  the 
circle  of  colatitude  7  is  given  by 
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2  * 

-4a     sin  VTr«sin  y^-^^ 

n 


_        .         ,,    ,     ^„„2        „,„2      „„„    /,,    ,.  \    27rT„,„2      „,„    ^,,    ,. ,    2Tr 
M- 


>        P^[-cos  -)/  -  sin  -y   cos    (k-j)  ^]sin  7    sin    (k-j)  — 


n  OP  m  Q  0^ 

-   H  IIIT""  C    (m,A)[P^(cos  7)]^msinin(k-j)  ^   . 

In  this  expression  each  of  the  sums  with  respect  to  j  has  the  form 

r^  ^[cos  (k-j)  ^]sin  (k-j)m^  . 

This  sum  is  equal  to  zero  and  so  we  find  that 

The  velocity  to  the  east  of  any  one  of  the  vortices  on  the 
circle  0  =  7  is  given  by 

4a  sin  vw^,     =   -4b  sin  vttP  (cos  7) 

+  n^[P^(-cos  -y)  -C^(0,A)P^(cos  7)] 

+  ^L  J~  P'[-cos  7  -  Sin  7   cos   j  —]2   cos  7    sin   j  - 

-   ^x  YIZ  HZ  ^v^"^^^^^    ^^°^  'y)F!!(cos  7)    cos  jm  rj 
j=l  m=0 

o 

=  4a     sin  V7r«0    . 
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We  have  now  shown  that  if  at  time  t  =  0  one  vortex  of 
strength  [i      is  at  the  north  pole  and  the  others  of  strength  p. 
occupy  the  positions        -  :    ■   ,  . 

i^y^j)    =  ((j-1)  ^'   7)  .    J  =  1,2,. ..,n  , 

then  the  polar  vortex  stays  at  the  north  pole  while  the  others 
remain  on  the  circle  of  colatitude  6   =  -y   and  each  moves  along  it 
with  constant  angular  velocity  O.   The  strengths  |j.^  and  \i   can  be 
chosen  so  that  0  =  0.   On  the  other  hand,  if  \i^   and  (i  are  pre- 
scribed it  may  be  possible  (depending  on  the  boundary  condition) 
to  choose  b  and  7  so  that  the  vortices  are  stationary. 
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5.  Linear  Stability  Equations 

The  motion  of  the  vortices  described  in  the  last  section  is 
given  by 


(5.1) 


(*,,0J  =  (^,,0) 


i^y^j)    =  ((J-1)  ^+Ot,7)  ,    j  =  1,2,. ..,n  . 


In  order  to  investigate  the  stability  of  this  motion  consider  a 
second  set  of  vortices  with  one  member  of  strength  \i      situated  at 

(5.2)  ('i'o'^o^  =  (a^+Ot,P^) 
and  with  the  others  of  strength  ^l  situated  at 

(5.3)  (*j'^j^  =  (a^. +  (j-l)  ^+ot,7+p^.)  ,    j  =  1,2,  ...,n  . 

Suppose  that  at  t  =  0  the  quantities  a.,  j  =  1,2,  ...,n;  and  p  , 

J  J 

j  =  0,1,2,  ...,n,  are  small.   If  these  quantities  remain  small  as 
time  goes  on  then  the  second  set  of  vortices  will  always  be  close 
to  the  first  set  (5.1).   A  condition  for  linear  stability  can  be 
fo\and  if  we  can  linearize  the  equations  which  result  when  the 
coordinates  of  the  second  set  of  vortices  are  substituted  in  the 
motion  equations  (4.10)  and  (4.11).   The  linearization  cannot  be 
performed  directly  because  a  need  not  be  small  in  order  to  have 
the  second  set  of  vortices  near  the  first  set.   However,  it  follows 
from  our  supposition  that 
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X  =  B   cos  a     a . 
o     o       o      J 


> 


>   j  =  1,  2,  .  .  .,n 


y  =  P   sin  a     P . 


are  small  at  t  =  0;  and  linearization  with  respect  to  these 
quantities  is  sufficient  for  the  purpose  of  developing  a  linear 
stability  criterion. 

After  somewhat  lengthy  computations  the  linear  stability 
equations  turn  out  to  be: 

(5.4)    qi(x^-Oy^)  =  (A^-B^^  sin^j  ^)y^ 


n 


2Tr  ,  ^' 


n 


2w 


+  C  (-y)  5   a.  cos  (j-1)  £1  +  C  (7)  ^        p.  sin  (j-l)  — 
.1=1  "^  .1=1   "^ 


n-1 


,2,  27r 


(5.5)    Qn  (y  +OX  )  =  -(A  -B  }         cos'^j  — )x 
\^   -jj        ^l^^o    o'     ^o   o  4— rr     "^     n   c 

J=0 


+  ^0(7)  XZ   «j  sin  (j-l)  ^-C^i-y)  \        ^.    cos  (j-l)  ^  , 


1^, 


&r 


27r- 


(5.6)   -q-L  sin  7.p^  =  ^  C^(7)[x^  sin  (k-l)  ^  -  Y^  cos  (k-l)  ^] 


n  n 

+  7Z  \.P.  +  ZII  Bj^Ja  -a  )  , 

JT^k 


j= 


(5.7)    q^  sin  7. a. 


_  1^  C'(7)[x  cos  (k-l)  — +  y^  sin  (k-l)  ^] 
(X   o^'^^'-o     ^    ^n   "^o     ^    'n 


n 


n 


>    A,  .a.  +  ^    C,  .P  .  +  Cp, 
^  "kj  J   4^  kj  J    ^k 

JA 


JA 


where 
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qi  = 


A    = 


4a      sm  VTT 


M- 


4b 


M-, 


^ilL^   P;(l)+-^   P;(1)[C^(0,A)-C^(1,A)P;(1)] 


+   n[C,(0,A)p'(l)P    (cos   -y) 


+   2C    (2,a)P"  (1)P"  (cos  7)  -  P^(-cos   7)    cos  7]    , 


B      =   sin^7[P"(-cos   7)  +4c^(2,A)P;;(l)P;;(cos  7)]    , 


0^(7)    =  P' (-cos  7)    sin  7  -  C^(l,A)P^(l)P^(cos  7)    sin  7   , 


\; 


27r 


p'(-cos  7,   .)    sin  7   cos  7   sin   (k-j)  — 

+   P"(-cos  7j^-)2   sin^7   cos  7    sin^(k-j)   ^  sin    (k-j)  — 

+   Sin  7  'YH  C    (m,A)Fj(cos  7)F^    (cos  7)111   sin   [m(k-j)  — ]    , 


m=l 


"kj 


.    2 


2Tr 


=  P^(-cos  7^j)    sin  7   cos    (k-j)  — 

„/  X       .    4         .    2,,      .X    27r 

+   p^(-cos  7      )    sm  7    sxn    (k-j)  — 

00  w.  00  Ptt 

-  }         C    (in,A)[F'  (cos   7)]      m      cos    [in(k-j)  —]    , 


m=l 


C,   .   =   -P 
kj  V 


(-cos  71^  •)•  [sin  7  +  cos  7   cos    (k-j)  — ] 


^kj 


+   P"(-cos  7j^.).sin   27  •  sin    (k-j)  ^ 

+   sinS  -IZ  C    (in,A)[i^'(cos  7)]^   cos    [in(k-j)  ^] 
m=0 
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=  ^^  ^in  v^  P"(cos  ^)  sinS 

+  -^   sinS[P;;(-cos  7)  +C^(0,A)P^(cos  -y)] 


+  ZZ   {Pv^-cos  7j)  cos  j  ^  +  P^(-cos  7_.)  sin^2-Y  •  sin^j  I] 
+  2~~  2ZI  C„(^'^)^  (^°^  -y)P^'(cos  7)  sinS  cos  mj  ^ 
+  ^  C^(m,A)[Pj;  (cos  7)]''  sin  7  , 


2   ,   .  2        /,   .  \  27r 
cos  7  .  =  cos  7  +  sm  7  •  cos  (k-jj  — —  ,   ■ 

2   ,   .  2        .  27r 
cos  7  .   =  cos  7  +  sm  7   cos  j  —  . 

The  linear  system  (5.4)- (5.7)  has  constant  coefficients. 

Therefore  the  investigation  of  the  stability  of  the  solution  of 

the  solution  of  the  system  can  be  started  by  substituting 

X  =ae^,a,  =a,e^,V  =be   ,6,  =b,  e    for  the  variables, 
o    okK'^o    o'k    k 

This  substitution  leads  to  the  characteristic  determinant  which 
theoretically  could  be  studied  in  the  usual  way;  but  it  is  clear 
that  the  required  analysis  for  the  general  case  would  be  too 
complicated  to  allow  the  deduction  of  a  practicable  formula  for 
linear  stability.   If  a  free  polar  vortex  is  involved  with  two  or 
more  other  vortices  it  is  best  to  subject  the  characteristic 
matrix  to  numerical  analysis.   C.  Leiva,  [2],  is  doing  this  and  he 
is  analyzing  the  results  for  the  case  in  which  the  cap  coincides 
with  the  northern  hemisphere. 

When  the  polar  vortex  is  either  absent  or  kept  fixed,  the 
system  (5.4)- (5.7)  simplifies  so  much  that  the  development  of  a 
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linear  stability  formula  becomes  feasible.   If  the  polar  vortex  is 
absent,  then  \i     =   0,      If  the  polar  vortex  is  present  but  con- 
strained to  remain  at  the  pole,  then  x  =  y  =0.   Hence  for  either 

^  o    o 

of  these  cases  our  system  reduces  to 

(5.8)  -qA  =  r:  \jpj  +  £:  \3^\ - °j)  > 

J   -^  <j   -^ 

(5.9)  ^A  =  -r:\f^^r:^Kf^^^\ 


where 


q      =    q,     sm   7    = 


4a     sin  vt  •  sin 


a 


2        ^1   ^^''    '  H 

For  reasons  stated  above,  let  us  devote  the  remaining  analysis  in 
this  paper  to  the  system  (5.8),  (5«9)  and  its  implications. 

Summations  show  that  (5.8),  (5 "9)  possess  two  invariants 
given  by 

n   . 


(5.10) 

and 

(5.11) 

where 

D 

=  ^  °-i 

0 


'2  5"^  =  '=^°' 5^^ 


^_l     r  -P^(-cos  7^.  )•  (sin  7  +COS  7   cos   j  -^) 
^^^    ^  +   P"(-cos  7.)    sin^27   sin^j  - 


+   sin  7  J        ^~  C    (m,  a)[F     (cos  7)]      cos  mj  —   . 
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From    (5,10)    and    (5. 11)   we    obtain 


n 


and 


n 


5 


11  \  =  "<J2«1 

k=l 


a^  =    (C+D)n^-^t  +ne. 


where  ^-,  and  |_  are  small  quantities.   If  we  set 


\   =   \^^2h 


we  see  that 
(5.12) 


"k  "  ^'k^  (C+D)^^t +^2 


n 


fci 


Pv  =  0  , 


(5.13) 


n 


£; 


\  =  0 


and  we  find  that  p,  and  a,  must  again  satisfy  the  same  equations 
as  p.  and  a  satisfy,  namely  the  system  (5.8),  (5«9).  In  matrix 
form,  we  have 


(5.14) 


+   q. 


^^ 


=  0  . 


Now  if  we  take 


r\j  ru 


\« 


st     ~    ~   st 

'    "k  =  \^ 


then  (5,14)  reduces  to 
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(5.15) 


V 


a     13\ 
C     a 


^ 


+  q2s 


=   0 


Here    the    subitiatrices   (X ,     f^    and     ^   are    right   circulant  nxn 
matrices.      It   can  be   shown   that   the    characteristic   roots   are   given 
t>y 

q.2^^i£)   =   -A(^)    +  /-B(^)E(i) 


(5.16) 


where 


i   =  1,  2,  . .  .,n  ; 


q„s„(^)    =   -A(^)    -   y-B{£}E{£) 


^2   2 


n  i(j-k) 

A(i)   =^\j- 
jVk 


ill 

n 


P^(-cos7    )sin7COS7Sinj  — 


n-l 
■i   ) /  +  P" 


•\ 


+  P"  (-cos  y  .)  '2  sin^-y  cos-i/  Sin  j—  sin  j  — 
V  '.  1  '  '  n  "^   n 


• 


CO 


+  )        C    (m,  A  )P^(cos  7  )P^    (C0S7)  sin^  .msinjin     ^ 
rn=l     V  V  V 


n 


sin   j&  ^ 


which  is  a  pure  imaginary;  while 

,27ri  1 


n 


n-l 


1  -  e 


£(j-k)- 


n 


1-.1  /  \    .    2  .  27r 

P   (-cos^.jsm  7C0S  J  — 


2Tr 


>         (1-cosj^^)    ;  +p;;(-cos7j)sin%sin^j  -^ 


CO 


2    2 


+  2 ^  C    (m,A)[F{|"(cos  7)]     m     cos  jm  -^ 


2jr_ 
n 


and 


31 


n  ^(j-k)^ 


E(^)    =   C   +  2^  C      e 
jVk 


n 


4b   sin  VTT   p..(cos^)slnS 

(X  V  '  ' 

+  _2  sinS[P;(-cos7)  +  C^(0,  A)P;;(cos'y)]  •       ■  ' 

^_1    /P^(-COS7    )[cosj^-    (sin^7+ COS^-yCOS  j^)cos  j^-^] 


^5 


+  P^( -COS  7j)  sin^27  sin  0^.(1+ COS  j£  -^)    ■ 

+  ^ZZ  m  C^("1'^)p!JJ'(cos  7)f)J    (cos  7)  sin  7  cos  jm-^ 
j=l  m=0 

+   sin  7  TZZl  C    ("^>^)[^    (cos7)]^cos  jm^  cos  ji  ^ 


are   real. 


Since   a,  ,    P,     are    supposed   to  be    small   at   t  =   0,    the   quanti- 


<^    oj 


ties  a,  ,  P^  must  be  small  at  t  =  0.   We  define  the  motion  given  by 
(5.1)  to  be  linearly  stable  if  a,  ,  '^,  remain  small  as  time 
increases.   Then,  in  order  to  have  linear  stability  it  is  necessary 
that  the  real  parts  of  the  roots  (5.16)  be  less  than  or  equal  to 

zero.   Since  A(i)  is  a  pure  imaginary,  it  follows  that 

I  - 
(5.17)  B(^)E(^)  >  0  ■     :-  ■ 

is  a  necessary  condition  for  linear  stability.   In  other  words, 
each  root  of  the  set  (5.16)  must  be  a  pure  imaginary  number.   With 
this  condition  satisfied  we  still  need  to  consider  the  possible 
appearance  of  multiple  roots  in  the  set  (5.16).   For  example  if 
^  =  n  then  A(n)  =  0;  B(n)  =  0;  and  consequently  there  is  at  least 
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a  double  root  present.   Presiomably  this  means  instability  since 
the  multiple  zero  root  indicates  solutions  of  (S.l-^)  which  have 
the  form  t"^"''".   However,  if  the  zero  root  is  just  a  double  root  it 
can  be  shown  that  solutions  a,  ,  P^^  which  in  part  depend  linearly 
on  t  are  precluded  by  the  invariants  (5.12)  and  (5.13). 
We  notice  of  course  that 

has  a  part  which  depends  linearly  on  t.   This  is  a  reflection  of 
the  fact  that  if  the  vortices  (5.I)  are  slightly  displaced,  then 
the  angular  velocity  0(7)  will  change  by  a  small  amount.   The 
quantity  (C+D)  is  equal  to 

C+D  =  4a   sin  vtt  •  sin  7  -v-  ■> 

Double  roots  other  than  a  double  zero  root  may  imply 
instability. 

If  [i,  =0  the  expectation  is  that  at  least  one  domain  of 

linear  stability  is  the  cap: 

0  <  7  <  A   <  A  . 
—  '     s 

In  fact,  if  li.  has  the  same  sign  as  |x,  and  the  polar  vortex  is 
kept  fixed,  then  an  investigation  of  the  signs  of 

B(i),  E(i),    £   =   1,  ...,n-l  ; 

shows  that  the  product  B(i)E(^)  is  positive  (no  matter  what  n  is) 
provided  7  is  sufficiently  small.   It  is  interesting  to  note  that 
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for  the  case  of  a  ring  of  Bessel  vortices  the  linear  stability 
condition  is  not  satisfied  for  the  neighborhood  of  the  origin  if 
n  >  8.   The  boundary  angle  0  =  A   can  be  approximately  determined 
by  computing  B{£)    and  E(^)  for  various  values  of  7.   For  any  given 
case  this  would  require  considerable  numerical  analysis.   As  an 
example  of  what  is  involved  we  turn  to  the  case  of  three  vortices 
on  a  circle  of  latitude. 
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6.  Three  Vortices  on  a  Circle  of  Latitude  in  a  Cap.   Polar  Vortex 
Fixed  or  Absent 

Our  intention  is  to  devote  the  rest  of  this  paper  to  a 
discussion  of  the  case  of  three  vortices  of  equal  strength  [i 
symmetrically  situated  on  a  circle  of  colatitude  0  j<  0  <  A.   When 
it  is  admitted,  a  polar  vortex  of  strength  |j,   is  assumed  to  be 
held  fixed  at  the  north  pole.   We  are  interested  in  the  linear 
stability  of  motion  of  the  vortices  whose  positions  at  any  time  t 
are  given  by 

(({>^,ej)  =  ((j-i)  ■^+ot,^)  ,   J  =  1,2,3  . 

Our  primary  object  here  is  not  to  present  an  extended  niomerical 
analysis;  but  rather  to  record  the  quantities  which  need  to  be 
evaluated  for  such  an  analysis. 

The  angular  velocity  of  the  vortices  in  the  equilibrium 
position  is  given  by 

(6.1)     ^a^sin  VTT  .^  ^  -4b    sin  vtt  .  p^(cos  ^) 

+  -^   [P^(-cos   -y)  -C^(0,A)P^(cos  7)] 


I  2         12 

+  3P„(-cos  7  +75-  sin  •y)«cos  -y 

C    (m,A)F^'(cos  7)I^(cos  -y)'(1  +  2   cos  m  -f) 


V 
00 

LI 

m=0 


For  the  present  case  of  three  vortices  we  find 
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B(l)    =   B(2)    =   B    , 
E(l)    =  E(2)    =   E 


where 


(6.2)      B  =  -  ^  P' (-cos^-y  +-i  sin^^)    sin^' 


2     V 


7 


"  Q  ?  1  ?  U- 

"•"  i  P^(-cos  -y  +^  Sin  ^)    sin  -y 


00 

-  5  YZ  C    (m,A)[Fj^(cos  7)]2  m^  cos  m  -^ 
rn=l     V  V  :? 


and 


//r  -J  \      -c.       4b   sin  VTT      t,,,  /  \      .2 

(60)      E  = •  P^(cos  7)    sm  7 

+    -J-    [P;;(-COS    -y)   +C^(0,A)P;;(COS    -y)]      Sin^^ 

3    Bi  /  2^1.2.  2 

-  ^   P^(-cos  "Y  +-^  sm  7)    cos  -y 

Q  21  2  2  '  ■ 

+  -^  P"  ( -  c  o  s  -Y  +  -^   s  in  -y  )    s  in   27    ■ 

+    sln^7  .y^  C    (in,A)l^(cos  7)P™    (cos   7)«(l  +  2   cos   m  ■^) 

p  OP  J  p  p 

+   sin  7  •  )         C    (m,A)[F     (cos  7)]    (l-cosm-^)    . 

For  £  =  '3 ,  B(3)  =  E(3)  =  0.  This  means  a  double  zero  root  in  the 
set  (5.16)  for  three  vortices.  There  are  no  other  multiple  roots 
and  hence   the   linear   stability   criterion   is 

(6.4)  B(^)E(^)    >   0   ,         £   =   1,2 


or 
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(6.5)  BE    >    0    . 

If    (6.5)   holds   then  the   component   periods   of   the  motion  about 
the   equilibrium  position   are   given  by 


s(^)    _  A(^)  ii/BE  .-12 

where 


(6.6)      iA(l)    =  I  P^(-cos^'Y  +i  sin^^)    sin  7   cos  7 

Q  2         12  3 

+  2-  P"(-cos  7  +-^  sin  7)    sin-^7   cos  7 


+  75  sin 7  '^     "  C^(m,A)F^(cos  7)1^    (cos  7)  msinm 


and 

A(2)    =   -A(l)    . 

Recall  that  if  the  boundary  condition  is 

^((t>,A)  =  c 

then 

, c 

^   -   P  (cos  Aj  ' 

V  ^ 

P  (-cos  a) 
Cv^°^^)  =  -  P'icos  A)  ' 

,  ^,  .   p"(-cos  A) 

c  (m,A)  =  -2(-i)"^ ;: " ::;!|  ~ >   ^^^ 


V 


If  the  boundary  condition  is 

then 
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c 

sin  A  P..  (cos  A) 


P'  (-cos  A) 


C  (0,A)  -  _r- —  > 

V         P^(cos  A) 

„  /     ,  -,  \   P^  (-cos  A ) 

If  the  cap  coincides  with  the  sphere  then  , : 

b  =  0  ;    C  (m,7r)  =  0  . 

V 

Some  ideas  about  the  effects  of  imposing  a  boundary  condition  for 
an  arbitrary  cap  can  be  drawn  from  a  comparison  of  the  motion  con- 
fined to  the  northern  hemisphere  with  the  free  motion  on  the  whole 
sphere.   If  the  cap  is  the  northern  hemisphere  then 


=  0=1' 


-  ol    ^\'^   r(v  -m-1) 
~  '^^~^'     r(v  +m+lj  ' 


e^  =  dK-y)     ^ — .   . .  \  ,    m  >  0 


and  the  series  which  appear  in  the  above  expressions  can  be  summed 


TT 


by  using  (4.7).   We  find  that  for  A  =  -^  : 
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(6.7) 


4a     sin  vtt     _        -^b 


M- 


O 


SXn     VTT  -r,! 


u- 


.  P^(cos   ^) 


M-, 


+  _o    [p'(-cos  7)  ±  P'(cos   7)] 
la         V  V 


P         1  2 

+  3P   (-COS  7  +-^  sin  7)    cos  7 


?         1  2 

±    [2P'  (cos   27)  +  P    (cos  7  +-^  sin  7)]    cos  7    , 


(6.8)      iA(l)   =  I  P^(-cos^7  +-i  sin  7)    sin  7   cos  7 


Q  2         12  3 

+  4-  P"  (-cos  7 +-5  sin  7)    sin^7   cos  7 


^  r^^-cub  7  -r-^ 


+    [P^(cos^7  +-^  sin  7) 

-  i  P"(cos^7+-i  sin^7)    sin^7]  4  sin  7  cos  7    , 


A(2)    =  -A(l)    , 


(6.9)      B 


I  P^(-cos^7  +^  sin^7)    sin  7 

+  2-  P"  (-cos  7  +-^  sin  7)    sin  7 


I  2  12 

±    [P   (cos  7  +-J  sin  7) 


I  P"(cos^7  +i  sin^7)    sin^7]   |  sin  7 
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(6.10) 


_        4b    sin  VTT   T^ii  /  n      .    2 

E  = P^(cos  7)    sm  7 


M- 


+  -^  [p;(-cos  7)  :pp;;(cos  7)]  sm^^ 


-  ^  P^(-cos  7+2   sin  7)    cos-^7 


Q  2         12  2 

+  rj?-  P"  (-cos  7+75  sin  7)    sin   27 


p'(cos   27)2   sin^7+P"(cos   27)2   sin^27'' 

■n'  /         2      ,  1       .    2    w„      1  2    ^         \ 

-   P^(cos  y  +2    211^  7)  (2  -"2   cos   7)        , 

+   P^(cos^7+i  sin^7)  .^  sin^27 

1 

If  the  boundary  condition  is 

the  upper  sign  is  to  be  taken  with 


If  the  boundary  condition  is 


fgi^  ,^)    =  o 


the  lower  sign  is  to  be  taken  with 


b  =  - 


pToT 

V 


A  numerical  analysis  for  the  caps 
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0     <    0     <    TT 

amd 

0  <  e  <  I 

will  be  presented   in   a  separate    report  by   C.    Leiva   [2] . 
For   the   arbitrary  cap 

0   <   0    <  A    <  I 

each  series  which  appears  in  the  expressions  (6.1),  (6.2)  and  (6.6) 
can  be  approximated  by  the  sum  of  its  first  n  terms.   Various 
methods  can  be  used  to  estimate  the  error  involved.   One  method  is 
to  use  the  formula 


-m/2      1 


dt 


X 


Py^^)        ,l-X,m/2 

fim+ry  *  ^  i+x ' 


x<t<l;        m^l 

and  the   formula 

(^    1o^  f    n  ^m  r(v-m+l)    r    ^\^u         ^v^^^       ,l-x^m/2         1 

^6-12)  (-1)    r(v+m+l)   ^-K   ^^^^  =^rWJ'  ^T^^       'Y:^' 

m  _>  2 
in  conjunction  with 

,^    >  ,   xm  r(v  -m+1)  _      cosh  qrr 

(6.13)         (-1)  r(v+»+i)  --,j^(o)^'(o)  ■ 

V      V 

For  example,  consider  the  series  which  appears  in  the  expression 
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for  O,  namely 

m  =  -  ^  C  (m,A)P^'(cos  7)r^(cos  7)(l  +  2  cos  m  ^) 
in=0  V       V  V  ;j 

Let  the  boundary  condition  be        .       ■       ■   ■      ' 

^((t),A)  =  c  . 


so  that 


P  (-cos  a) 
V  ^ ' 

P  (cos  Aj 

V 


C  (m,A) 


V 


If  Sp,  the  sum  of  the  first  two  terms,  is  used  to  approximate  ^ 
we  have 

f-    _  ,   ^     OO  P5'"(-C0S  A)  P^'^'CCOS  7)P5'"(C0S  -y) 

D  cosh  qir  -r y  ^ v       '  ^  v       " 


2"^-       ^     fc2p5"^(cosA)       P5"^'(0)p5"^(0) 


Then  from   the   formulas    (6,11)    and    (6,12)   we   find 


1  '       ^  +2   A    ,      2  7,3ra 

P    (cos    A)[P    (0)]"^  m=2  ^  ^ 


Tr(S„  -  >        )  o      00  o  -x  c        3m 

>^~~        <   P    (-cos   A)[P   (cos  -y)]^  T~   (cot^  4  tan^  ?) 
6   cosh   qir  v  ' '-    v  ^r^  2         ,       2 


Since 


^   =   cot  A  tan  ^  <   1 


we   obtain 
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c 


12  ^     -jriS^-YZL)         P^(-cos  A)[P^(cos  7)]2;;12 


T^  /     ^\rT3  /■^m2  ^^^  To  5  cosh  q-fT  t  To 

P^(cos  A)[P^(0)J    1-C  ^  1-C 


This  shows  that  ^   -^  oo  as  C  -*  1j  which  means  that  O  is  large 
when  the  vortices  on  the  latitude  circle  are  close  to  the  boundary 
circle.   It  also  shows  that 


as  7  -*  Oo   A  similar  analysis  can  be  applied  to  the  infinite  sums 
in  (6.2),  (6.3)  and  (6.6). 


43 


'         References      ,   i  ■,  , 

[1]   Morikawa,  G.  K.,  Geostrophic  Vortex  Motion,  Journal  of 
Meteorology,  Vol.  17,  No.  2,  April  I960,  pp.  148-158. 

[2]   Leiva,  C,  Courant  Institute  of  Mathematical  Sciences, 
New  York  University.   Report  to  be  published.  ■'-  - 

[3]   Morikawa,  G.  K.,  and  Swenson,  E.  V.,  Interacting  Motion 
of  Rectilinear  Geostrophic  Vortices,  the  Physics  of 
Fluids,  Vol.  l4.  No.  6,  June  1971,  PP.  1058-1073. 

[4]   Peters,  A.  S.,  Motion  of  a  Vortex  in  a  Rotating  Layer 
of  Fluid  Between  Rigid  Concentric  Spheres,  Courant 
Institute  of  Mathematical  Sciences,  New  York  University, 
IMM  398,  September  1973* 


44 


SECURITY   CLASSIFICATION   OF   THIS  PAGE  (When  Data  Entered) 

REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

1.     REPORT  NUMBER 

IMM  4o4 

2.   GOVT   ACCESSION   NO. 

3.     RECIPIENT'S  CATALOG   NUMBER 

4.     TITLE  (and  Subtitle) 

GEOSTROPHIC   VORTICES    ON  A   CIRCT,E    OF 
LATITUDE    IN  A   CAP   ON  A   ROTATING    SPHERE 

5.     TYPE  OF   REPORT  &  PERIOD  COVERED 

Technical 

6.     PERFORMING  ORG.    REPORT  NUMBER 

7.     AUTHORfsJ 

A.    S.    Peters 

8.     CONTRACT  OR  GRANT  NUMBERfs; 

AFOSR  7^-2728 

9.     PERFORMING  ORGANIZATION   NAME  AND  ADDRESS 

Courant    Institute   of  Mathematical   Sciences 

New  York  University 

251   Mercer   St.,    New  York,    N.Y.    10012 

10.     PROGRAM   ELEMENT.  PROJECT,   TASK 
AREA  a   WORK   UNIT   NUMBERS 

ARPA   Order   2774 

11.     CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force   Office    of   Scientific   Research 
l400  Wilson  Boulevard 
Arlington,    Virginia      22209 

12.     REPORT  DATE 

July   1974 

13.     NUMBER  OF  PAGES 

44 

14.     MONITORING  AGENCY  NAME  »   ADORESSf// d///eren(  Irom  ControtUng  Office) 

Advanced  Research   Projects   Agency 
l400  Wilson  Boulevard 
Arlington,    Virginia      22209 

15.     SECURITY  CLASS,  (of  this  report) 

not   classified 

15a.     DECLASSIFI  CATION/ DOWN  GRADING 

^"""^^^            none 

16.     DISTRIBUTION  STATEMENT  (of  this  Report) 

Distribution   of   this   document    is   unlimited. 

17.     DISTRIBUTION  STATEMENT  (ol  the  abstract  entered  In  Block  20,  It  dlllerent  from  Report) 

none 

le.     SUPPLEMENTARY  NOTES 

none 

t9.     KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 

none 

20.      ABSTRACT  (Continue  on  reverse  side  If  necessary  and  identity  by  block  number) 

This    report   contains   an   analysis   of   the   linear   stability   of 
the  motion   of  n+1  geostrophic   vortices   confined   in   a  polar   cap   on 
a  rotating   sphere.      It   is    supposed   that   initially   one   vortex   is   at 
the  north  pole  while    the   others    are    symmetrically   arranged   on  a 
circle   of   latitude;    and   that   the  motion   is   subject   to   a  certain 

DD     1    JAN    73     1473  EDITION   OF    1    NOV  65  IS  OBSOLETE 

45 


SECURITY   CLASSIFICATION   OF   THIS   PAGE  flWien  Data  Entered) 


SECURITY  CLASSIFICATION  Of  THIS  PACE(W7l«o  Data  Bnfnd) 


20.  Abstract 


condition  along  the  latitude  circle  which  bounds  the  cap.   The 
development  stems  from  the  demonstration  that  a  concentrated 
geostrophic  vortex  in  a  rotating  spherical  layer  is  characterized 
by  a  singular  spherical  harmonic  of  degree  v  and  order  zero. 
This  function  has  a  role  analogous  to  that  of  the  Bessel  function 
which  characterizes  a  rectilinear  geostrophic  vortex  in  a 
rotating  planar  layer.    


This  book  may  be  kept 


;j:?4 


FOURTEEN    DAYS 


A  fine  will  be  charHcd  for  each 

day  the  book  is  k 

ept  overtime.          \ 

A     -Cl 

^ 

■,0\^ 

GAYLORD    I<2 

PRIKTED    (N    U    S    * 

46 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfHTien  Dmtm  Bntend) 


Peters 


c.l 


Geostrophic  vortices  on  a 
circle . . . 


IMM-404 
Peters 


C.l 


AUTHOR 


GenstrnpTnic    vortices    on    a 


N.Y.U.  Courant  Institute  of 
Mathematical  Sciences 

251  Mercer  St. 
New  York,  N.  Y.    10012 


